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We present new upper critical field Hc2{T) data in a broad temperature region 0.3K < T < Tc 
for LuNi2B2C and YNi2B2C single crystals with well characterized low impurity scattering rates. 
The absolute values for all T, in particular _ffc2(0), and the sizeable positive curvature (PC) of 
Hc2{T) at high and intermediate T are explained quantitatively within an effective two-band model. 
The failure of the isotropic single band approach is discussed in detail. Supported by de Haas van 
Alphen data, the superconductivity reveals direct insight into details of the electronic structure. 
The observed maximal PC near Tc gives strong evidence for clean limit type II superconductors. 
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The discovery |]l],|j of superconductivity in transition 
metal borocarbides has generated large interest due to 
their relatively high transition temperatures Tc '^ 15 to 
23 K and due to the relation between the mechanisms 
of superconductivity in these compounds, in cuprates, 
and in ordinary transition metals. Another highlight is 
the coexistence of magnetism and superconductivity in 
some of these compounds containing rare earth elements 
P-p|. A study of the non-magnetic compounds such as 
LNi2B2C, with L=Lu,Y,Th,Sc p, is a prerequisite for 
the understanding of their magnetic counterparts. Ex- 
perimental data for LuNi2B2C Q] demonstrate beside 
a maximal positive curvature (PC) of Hc2{T) near T^ 
observed also for YNi2B2C l^,^,^, a weak T-dependent 
anisotropy within the tetragonal basal plane and a T- 
indcpendent out-of-plane anisotropy of the upper criti- 
cal field Hc2- Both anisotropics have been described |^] 
in terms of nonlocal corrections to the Ginzburg-Landau 
(GL) equations. In this picture the PC of H^2 {H\\ to 
the tetragonal c-axis) is caused, almost purely, by the 
basal plane anisotropy. However, it should be noted that 
the reported anisotropy of Hc2 for YNi2B2C is signif- 
icantly smaller than for LuNi2B2C |l7|,p|JlO| whereas its 
PC is comparable or even larger. Further explanations of 
the unusual PC of Hc2{T), such as quasi-2D fluctuations 
PI, are excluded by the underestimation of Hc2{T) at 
low-T [^ and the observed weak anisotropy. The quan- 
tum critical point scenario |13| as well as the bipolaronic 
one |l^ can be disregarded because the slope of Hc2 {T) 
decreases for T— >0 (see Fig. 1). Local density approxi- 
mation (LDA) band structure calculations ^4 15 1 predict 
a nearly isotropic electronic structure with rather com- 
plicated bands near the Fermi level Ep. However, in 
analyzing the superconductivity in terms of an isotropic 
sm^Ze-band (ISB) Eliashberg model, the multi-band char- 
acter and the anisotropic Fermi surface have been widely 
ignored so far. 

Here we present and analyze theoretically new data of 
Hc2{T) in a broad interval 0.3K< T <Tc for high purity 



LuNi2B2C and YNi2B2C single crystals. We show that 
typical features of both compounds, such as i/c2(0) ~ 

8 to 10 T and the unusual PC of Hc2{T) for T - 0.5Tc, 
cannot in any way be explained consistently with the nor- 
mal state properties within the ISB approach. Instead we 
propose a two-band model (TBM) approach. To clarify 
its relationship to the extended saddle-point model |19] 
which also predicts a PC is beyond the scope of our letter. 
Platelet shaped LuNi2B2C and YNi2B2C single crys- 
tals with a mass of ^ 1 mg were grown by a high- 
temperature flux technique with Ni2B as flux. The values 
of Tc, 16.5 K and 15.7 K, have been determined by low- 
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FIG. 1. Experimental data for Hc2{T) of LuNi2B2C (mag- 
netic field H\\ the c-axis) compared with theoretical curves; 
(i) the isotropic single band (ISB) model with Fermi velocity 
uf=2.76T0 cm/s and various impurity scattering rates -yimp 
explained in the text and the legends and (ii) the two-band 
model (TBM) with «f« (i=l,2) in units of 10^ cm/s. 



field ac susceptibility xC^) with transition widths AT^ 
= 0.2 K. The upper critical field 7?c2(T) along the c- 
axis, shown in Figs. 1 and 2, has been measured re- 
sistively for fixed T adopting the midpoint criterion: 
p(7?c2,T)=0.5p(i7==0,r=17K)= 0.5p„. The transition 
width Ai?=[i/(/9=0.9p„) — H{p=Q.lpn)] increases up to 
0.75 T (1.5 T for LuNiaBaC) at T < 0.5 K start- 
ing from a nearly constant value of 0.3 T at T >4 
K (8 K for LuNiaBaC). The low residual resistivity 
p(0) « p„=2.5^rJcm and the ratio p(300K)/p„ =43 (27 
for LuNiaBaC), together with the observations of mag- 
netoquantum oscillations |T^-Q9] indicate a high quality 
and a low impurity content of our samples. This suggests 
that we are in the clean limit in terms of the traditional 
theory of type II superconductors (20| . In this limit one 
has to consider the electronic structure in more detail. 
We restrict ourselves to an effective two-band model ||2l| ] 
which, especially on the simple BCS-level, has a long his- 
tory [g2|. Due to the neglect of strong coupling effects, 
a BCS-like theory is not expected to describe real super- 
conductors quantitatively. Such effects must be studied 
within the Eliashberg theory [pSJ-pTJ. To calculate i?c2 
we have solved numerically the corresponding linearized 
equations of Ref. 23, 

Cbi in) = LOn + ttT > [\.j {m — n) + 

j,'m 
+Smnhimp-i,3 + ■y-rap:,j) /2TTT]sgn{uJm) , (1) 

A,j(n) = 7rr^[A.ij(m - n) - p*5ij0{ujc- \i^m\) + 

j,m 
+Smnhtmp;Lj - lLp;t,j) /'^T^T]Xj{m)Aj{m) , (2) 



X{n) = (2/VA) / dqcxp{-q^) 
Jo 

X ta.n^^ {qy^i/ {\u}.,{n)\ + i/ZB-ffc2Sgn(a;„)), (3) 

A = eH,2vlJ2, (4) 



A,,,(n) = / dujujal^F{uj)/{uj^+ujl). 



(5) 



The bands at Ep are labeled by i,j. Here w„=27rr(2n-|-l) 
are the Matsubara frequencies, afF{uj) and A^, denote 
the spectral density, the superconducting order param- 
eter of the j*'* band respectively. In our approach, as 
in any two-band model, two gaps, below and above the 
BCS-value of 3.5kBTc, occur naturally. In general, in- 
terband coupling (i^j) mediated by phonons (a? -F(aj)) 
and impurities is important. Since there is no experi- 
mental evidence [^,^ for the presence of magnetic im- 
purities in high quality samples, we neglect the magnetic 
scattering rate 7|„p. For the quantification of the non- 
magnetic counterpart ^imp ~ 27rTD, the Dingle tem- 
peratures, Tjj, measured by the de Haas- van Alphen 
(dHvA) effect are very suitable (l7-19|. The experimen- 
tal values Td=2.8 K and 4K reveal 7imp=18K and 25K 
for our YNiaBaC and LuNiaBaC single crystals, respec- 
tively, indicating that the clean limit is reached since 



limp f^ 2Ao ~ 51K holds for both samples, where 2Ao 
denotes the smaller of the two gaps. Hence, the scat- 
tering by impurities can be neglected setting jimp=0- 
In the weak coupling limit of an ISB case, Eqs. (1-5) 
are equivalent to the well-known WHH-theory pQ]. Any 
anisotropy of Hc2 can be described by a similar, but much 
more tedious, system of equations |26|. Since the mea- 
sured anisotropy is relatively weak, it will be ignored for 
the sake of simplicity. Therefore, only H^2{'^y^ will be 
compared with those computed for our isotropic models. 

The standard ISB model | p7| describes quantitatively 
the renormalization of the physical properties of metals 
due to electron-phonon (el-ph) interaction. The input pa- 
rameters of the ISB model are the density of states at Ep, 
N(0), the Fermi velocity vp, the impurity scattering rate 
7imp, the Coulomb pseudopotential /i*, and the spectral 
function a'^F{uj) of the el-ph interaction. These quan- 
tities can be determined from a few experimental data: 
the normal state low-T electronic specific heat 7sT, the 
plasma frequency ujpi inferred from the optical conduc- 
tivity, iJt,2(0), Tc and its isotope exponent a, as well as 
the normal state low-T dc resistivity p(0) w p{Tc) which 
similarly to Tp, gives a direct measure of the sample pu- 
rity. We adopt for the Coulomb pseudopotential ^*=0.1 
and hoJc=QOO meV for the energy cutoff in Eq. (2). The 
total el-ph coupling constant, A=2j du)a'^F{io)/u), can be 
estimated from the boron isotope effect as ~0.2 [pi| and 
the well- resolved phonon spectrum |3^ . 

We first consider LuNiaBaC. To find a lower bound for 
A, we accounted for only the high-energy carbon phonons 
centered at 50 meV and the boron branch at 100 meV. 
Fitting the experimental ub and Tc values, we obtained 
the partial coupling constants Aioo=0.31, A5o=0.22, and 
A=Aioo + A5o=0.53, where the subscripts denote the cor- 
responding phonon energies in meV. An upper bound of 
A=0.77 has been found using the Lu phonons centered 
near 9 meV (A9=0.34) and the same B band (Aioo=0.43) 
as in the case before. In the following a wide averaged 
spectrum with A=0.65 (Aioo=0.37, A5o=0.12, Ag^O.lB) 
will be used which reproduces the experimental values 
of ub and Tc- A^(0)=11.8 mJ/mol fc^K^ has been esti- 
mated from the experimental value of 7s = 27r^A;g(l-|- 
A)7V(0)/3= 19.5 mJ/mol K^. The value of wf=2. 76-10^ 
cm/s follows from the experimental value |33| of the 



plasma frequency Tujjpi = \J ^Tre^v^pN {{)) /i =4.0 eV. The 
analogous values for YNizBaC are A=0.637, iV(0)=ll.l 
mJ/mol fc|K2, vp^ 2,-W cm/s and i?c2(0)=2T, where 
the data of Refs. 8 and 34 have been used. 

We solved Eqs. (1-5) with these parameter sets for 
two types of spectral densities a^F(cj): (i) a wide aver- 
aged spectrum and (ii) a single Einstein mode peaked at 
huJE=4:2A meV chosen to yield the experimental Tc=16.5 
K for LuNiaBaC using the same value of A=0.65 as in the 
first case. The results are shown in Fig.ll|. Note that in 
the intermediate coupling regime under consideration, as 
expected, Hc2{T) is insensitive to details of the shape 
of a^F{uj) p^] and, in the clean-limit case, also insen- 



sitive to the actual value of the small scattering rates. 
Comparing the LuNi2B2C data with the ISB curves one 
clearly realizes strong deviations. In particular, there is 
a discrepancy of about 3 between experimental and ISB 
model values of ifc2(0). For YNi2B2C the discrepancy 
reaches even a factor of 5. Therefore it makes no sense 
to discuss any details of the shape of Hc2 (T) , such as the 
PC, resulting in deviations of the Hc2{T) curves of the 
order 10 to 20 % , until the reason for the large failure 
to account for the magnitude of i?c2(0) has been eluci- 
dated. We remind the reader how this serious difhculty 
was circumvented in previous studies. At first, frequently 
the quantity hc2iT) = -Hc2{T)/[T,{dH,2/dT)T=T^], de- 
scribing the shape of the Hc2{T) curve, has been con- 
sidered, but the absolute values of Hc2{T) have not been 
discussed at all m. At second, since within the ISB- 
model i?c2(0) is a monotonically increasing function of 
the impurity content, in principle, large iJc2(p)-values 
might be obtained. To check this approach Q we cal- 
culated the impurity scattering rate ^imp which is re- 
quired to increase Hc2{0) up to the LuNi2B2C value of 
7.6 T. Thus we obtain ^mp w 427 K which would lead 
to p(0) «17 /iJ7cm which strongly deviates from the ex- 
perimental value p{0) w2.5 /xJ7cm. In this context we 
note that our data and those of Ref. 35 show dependen- 
cies just opposite to those predicted by the ISB-model: 
Hc2{0) and Tc increase when ^imp decreases! At third, 
a further empirical parameter, the clean limit coherence 
length, has been introduced in Ref. 35. However, this 
results in the overdetermination of the model parameter 
set and the consistency of the two values vp obtained 
using (i) a clean-limit coherence length and (ii) normal 
state data has not been checked. Thus the ISB approach 
fails to explain simultaneously the three values 7s=19.5 
mJ/mol K^, hijpi^A eV, -ffc2(0)=7.6 T in the clean hmit, 
and the four values 75, ojpi, iJc2(0), p=2.5 /iOcm in the 
dirty limit. In addition, the ISB model is also unable to 
explain the small gap values 2AQ/kBTc <3.5 observed in 
microwave [Q, tunneling |10|, and dHvA [Q measure- 
ments. Furthermore the PC of Hc2{T) near Tc and the 
extended quasi- linear behavior of Hc2{T) down to T ~ 
1 to 2 K cannot be described within the ISB. The ISB 
model with a single vp contradicts the dHvA data which 
clearly show the presence of at least six different sections 
Fa, ...Fjj [Q|9| with roughly two or three groups of vp's. 
Turning to our TBM we solved Eqs. (1-5). For 
the sake of simplicity, the same phonon spec- 
trum as in the ISB case discussed above has 
been adopted. We achieve an excellent agreement 
with the LuNi2B2C data (see Fig. 1) for Ai=0.51, 
A2=A2i=0.2, Ai2=0.4, ^l^fj,*=0.1, WF2=3.7-10^cm/s, 
WF,i=0-96-10'^cm/s. For YNi2B2C we used the follow- 
ing set: Ai=0.5, A2=A2i=0.2, Ai2=0.4, ^|=/i5=0.1, 
WF2=3.8-10'^cm/s, Ui.i=0.85-10'^cm/s, A^(0) = 11 mJ/mol 
k^K^. Reproducing Tc=15.6 K, the adopted values of 
iV(0) agrees well with the LDA-value |1^ of 9.5 mJ/mol 
k^K^. The plasma frequency hujpi =4.4 eV is in accord 
with ?ia;„;=4.25 eV obtained in Ref. 34. From the ob- 



tained vp's it is concluded that Upi is mainly related to 
the second weakly coupled band. Then transport, opti- 
cal and tunneling data mainly exhibit the properties of 
that band whereas the strongly coupled band remains 
almost hidden. The calculated value of the penetration 
depth at 4.2 K, 100 nm, is in agreement with the data of 
Ref. 35. Our 7^=17.2 mJ/mol K^ should be compared 
with 75=18.2 mJ/mol K^ reported in Ref. 8. Finally, 
we arrive at 7Jc2(0)=9.4 to 9.9 T in good agreement 
with our experimental value -ffc2(0)=10.6 ±0.2 T. The 
experimental Hc2{T) curve of YNi2B2C together with re- 
sults of the TBM are shown in Fig. 2. The values vp ~ 
4.2-lO^cm/s and 0.7 to 1.3-lO^cm/s of the extremal orbits 
F^ and FVl/2, respectively, derived from earlier dHvA 
data [|l8|,^ on the same YNi2B2C single crystal do not 
deviate much from the parameters vp2 ~3.8 -lO^cm/s 
and vpi « 0.8 to 0.96T0^cm/s introduced empirically in 
our approach. The remaining deviations might be due 
to natural differences between vp^s on extremal orbits 
seen in the dHvA-experiments and the corresponding ef- 
fective quantities of our TBM which contains implicitly 
information on the whole Fermi surface. 

Further calculations within our TBM reveal that the 
PC of Hc2{T) mainly depends upon the strength of the 
interband coupling (A12, A21) and to a lesser extent upon 
the ratio of the Fermi velocities and the intraband cou- 
pling strength (Ai, A2). These findings may also be inter- 
preted in terms of the flat Ni-derived band near ¥ip and 
the dispersive bands crossing E^? seen in the LDA-band 
structure ||l4| . The strong mixing of these bands might 
be viewed as the microscopic origin for significant inter- 
band coupling. Further work is required to clarify this 
point. For low T <4 K, Hc2{T) is very sensitive to the 
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FIG. 2. Temperature dependence of Hc2{T)^ for YNi2B2C. 
Experimental points for the magnetic field H\\ c-axis. The 
Fermi-velocities of the two-band model (TBM) are given in 
units of 10^ cm/s. For the interaction constants see the text. 



actual vf,2/vf,i ratio. This is illustrated by the addi- 
tional curves shown in Fig. 2. The variation of vpi and 
vf2 results in slightly higher and lower Hc2 (0) values, re- 
spectively. Our empirical parameter sets for LuNi2B2C 
and YNi2B2C differ almost only in their vpi values which 
roughly scale with the Ni-Ni distance as d~^\_j^^. In this 
context the study of Sc(Th)Ni2B2C crystals having much 
reduced (increased) Ni-Ni distances is of interest. 

Our analysis reveals that both compounds can be well 
described within an effective two-band model provided 
there are at least two groups of electrons having (i) sig- 
nificantly different Fermi- velocities, (ii) strong coupling 
in the small-VF band, as well as (iii) sizeable coupling 
between the small-WF and the large-Wi? band. This case 
differs from the situation considered in Refs. 24 and 26. 
There, the stronger coupling is in the large-w^? band and 
the curvature of Hc2{T) near Tc is negative. A PC would 
only appear at intermediate T if the interband coupling 
and the impurity scattering are both weak. In contrast, 
in this region Hc2{T) shows almost no curvature in our 
model. In other words, the result of Refs. 24 and 26 can 
be understood as an average over two weakly coupled su- 
perconductors, the first with a high iJc2(0) but a low Tc 
and the second one with a small iJc2(0) but high Tc. In 
our case the isolated small- ii^ subsystem would have high 
values of A^iJc2(0), and Tc. The values of Hc2{0) and Tc 
of the coupled system are reduced by the second large- 
vp subsystem with weak interaction parameters. In this 
case, which to the best of our knowledge has not been 
considered so far, the PC of the resulting Hc2{T) near 
Tc becomes a direct manifestation of that interband cou- 
pling. In our TBM, the PC of Hc2{T), as well as Fc2(0), 
and Tc are suppressed by growing impurity content and 
the PC vanishes upon reaching the dirty limit with Tc ~ 
II K. Thus the observed maximal PC of Hc2{T) near Tc 
of multi-band systems of the type under consideration 
can be regarded as a direct measure of the high sample 
quality, as opposed to the widely spread belief attributing 
it simply to sample inhomogeneities. The latter scenario 
is excluded by the sharp transition in the x{T) data. 

To summarize, we have shown that the supercon- 
ductivity of pure non-magnetic borocarbidcs should be 
described within a multi-band picture. The isotropic 
TBM gives a reasonable starting point toward the under- 
standing of the mechanism of superconductivity in these 
compounds. Combined studies of quantum oscillations 
(dHvA) and Hc2{T) unified by Eliashberg analysis are 
found to be valuable supplementary tools to elucidate 
the specific role of subgroups of electrons having small 
wf's and strong coupling to bosons. These electrons may 
be readily overlooked by other experimental techniques. 
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